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Abstract 



On the simple model of interacting massless and heavy scalar fields it is demon- 
strated that the technique of heavy baryon chiral perturbation theory reproduces 
the results of relativistic theory. Explicit calculations are performed for diagrams 
including two-loops. 
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I. INTRODUCTION 



Heavy Baryon Chiral Perturbation Theory (HBCHPT), suggested in []I| is an important 
and effective method of calculation of different processes involving electro-magnetic and strong 
interactions (For review and references see 0). The authors of Ref. used the ideas of the 
heavy quark effective field theory which allowed them to avoid severe complications appearing 
in problem of relativistic treatment of baryons at low energies, encountered in . Jenkins and 
Manohar suggested to take extremally non-relativistic limit of the fully relativistic theory and 
expand in inverse powers of baryon mass M. 

In the heavy baryon approach one integrates out heavy degrees of freedom and expands 
resulting non-local operators in inverse powers of large mass. In terms of the relativistic per- 
turbation theory of the original field theoretical model (Feynman diagrams) heavy baryon 
approach corresponds to the expansion of integrands in the loop integrals in powers of 1/M 
with subsequent term by term integration of the resulting series |]^. The non-commutativity of 
the integration over loop momenta and the expansion in 1/M generates a problem of matching 
of heavy baryon approach to the original relativistic theory According Lepage's argument 
from the uncertainty principle, one should be able to compensate the difference between the 
results of "naive" heavy baryon and relativistic approaches by including additional terms into 
the Lagrangian of the heavy baryon approach. While the problem of this matching has been 
analysed at one loop level @J3fl' hest of our knowledge the matching procedure for 

higher order loops has not been studied. 

In the present paper we consider the matching problem on a two loop level on the example 
of the forward scattering amplitude in a scalar theory. The consideration of the non-zero spin 
and the non-zero transferred momentum makes calculations more tedious and less transparent, 
bringing nothing new and essential in the problem considered. In our calculations we use the 
technique of calculation of loop integrals by dimensional counting, developed in |[T^ . 

We explicitly show that heavy baryon approach reproduces the results of the original rela- 
tivistic theory at two-loop level. 



II. ONE-LOOP ANALYSIS 

Let us consider a field theoretical model described by the Lagrangian: 

C = _ 1m2$*$ - ^(j)df,d''(j) - g^*^(j) + Li (1) 

where $ is a complex scalar field with mass M, is a neutral massless scalar field, g is a. 
coupling constant and Li contains all counter-terms which are necessary to remove divergences 
(one can include also interactions with the derivatives and/or a larger number of fields and 
corresponding counter-terms) . 

To avoid complications due to the infrared singularities we work in six dimensional space- 
time. We use dimensional regularization and is a dimension of space-time. 

Heavy baryon approach to the processes which involve one heavy particle uses following 
expansion of the heavy scalar propagator (p^ = Mv^ + fc^, v"^ = 1): 
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FIG. 1. One-loop corrections to the heavy scalar propagator in the heavy scalar-massless scalar 
scattering process. 
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This expansion corresponds to the following Lagrangian: 

£i = -Mr ■ ^ + 1^) ^ (3) 

where the second term of £i is treated perturbatively. This Lagrangian can be obtained from 
the free part of C corresponding to heavy scalar field, defining $ = exp{iMv ■ x}ip. 

Let us start with the one loop self-energy correction to the scattering process in the original 
relativistic theory, depicted in FIG.l a) (The solid line corresponds to the heavy scalar and 
dashed line corresponds to massless scalar). The expression for this diagram is proportional to 
the following integral: 

f dJ^Q 
"^'^ " J [g2 + le] [{p + g)2 _ M2 + le] 

The straightforward integration yields {p' = Mv, p = Mv + l, v"^ = 1, l"^ = 0, p"^ = M^ + 2Mf ■/): 
+ r„-2) ^FAl.i'n/2;i-n;-—) (5) 



where pFg (oi, a^; 61, bq] z) are the (generalised) Hypergeometric functions of z |TT1 . 

Heavy baryon expansion for Jn is realized by expanding (according to Eq. (^) the integrand 
in 1/M and integrating the resulting series term-by-term. 

As it was observed in |T^, expanding integrand in powers of some parameter and changing 
the order of integration and summation one recovers that part of the value of the integral which 
can be expanded in powers of given parameter with non-zero coefficients. 

From the expression (^ we see that the first term can be expanded in powers of 1/M with 
non-zero coefficients, while the second one can not - it contains (M ) ^ . Hence we expect 
that heavy baryon approach reproduces the first term of the expression (Rl). 
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FIG. 2. One-loop corrections to the vertex in the heavy scalar-massless scalar scattering process. 
Indeed: 



J ^ [ <^ f 

J \Q^ + ie] \(v + a)^ - + le] J 



[g2 + ie] [{p + qf - + ie] J [q^ + ie] [{I + qf + 2Mv ■(/ + ?) + ^e] 
1 /• 



2MJ [q^ + re] . (/ + g) + + 
expanding the integrand and changing the order of integration and summation we obtain: 

J _ ^ { f d^'q 1 /■ d^'qiq + iy 1 

''^'^ ~2M\J [g2 + te] [v.{l + q)+ te] ~ 2M J [q2 + ^e] [v ■ {I + q) + tef J ^ " ' 

= - l) r (3 - n) i-2v . IT--' + - l) r (3 - n) {n - 2)i-2v ■ Z)-^ + 



(6) 



As it was expected, Eq. (|]) reproduces the expansion of the first term of the expression (|^). 

The second term of the Eq. which can not be expanded in powers of 1/M is analytic in 
momentum / and hence can be reproduced by introducing additional terms into the Lagrangian 
of the heavy baryon approach. 

Free propagators of the heavy scalar particle appearing in the expression for the considered 
diagram are apparently reproduced by heavy baryon approach. The same is true for all diagrams 
and we will not include the contributions of the free propagators in our analysis below. 

FIG.l b) schematically represents the first term in right hand side of the Eq. (or Eq.(^) 
FIG.l c) corresponds to the contributions of compensating terms. 

Next one-loop diagram we are considering here is drawn in FIG. 2 a) (one-loop vortex cor- 
rection to the light scalar-heavy scalar vertex). The result of this diagram is proportional to 
the following integral: 



J 
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d^'q 



[g2 + ie] [g2 + 2p'q + ie] [{p + q^ - + ie] 



-i(M^)' 7r2^ ; \ — ^3^2 1,1,3 -n/2; 2,5- n; ^ 



r (n - 3) - - ^ ' ' / ' ' ' J 

I 2\f-3 n /M2-p2\""^r(4-n)r(n/2-l) / , M2-p2^ 
-H^) i-3 ^3fJ./2-1,.-3;.-2;^ 
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On the other hand heavy baryon approach leads to: 

1 f d^'q 



J 



12HB 



4M2 J [g2 + ie] . (/ + g) + ie] [t, . g + ie] 
1 f d^'qiq + lf 



+ 



8M3 J [g2 + . (/ + g) + ^e]2 . ^ ^ 

The comparison of Eq. (|^) and Eq. (H) shows that heavy baryon approach reproduces that 
part of relativistic answer which can be expanded in inverse powers of M. The second part is 
analytic in / and can be reproduced by adding appropriate terms into the Lagrangian of the 
heavy baryon approach. FIG. 2 b) and FIG. 2 c) correspond to Eq. and the contributions of 
compensating terms respectively. 

The analysis of the rest of one- loop diagrams lead to the same result: heavy baryon approach 
reproduces those parts of diagrams which are non-analytic in the momenta and the remaining 
parts, analytic in momenta can be reproduced by adding terms into the effective Lagrangian 
of the heavy baryon approach. 

III. TWO-LOOP ANALYSIS 

Two-loop diagrams have more complicated structure. Let us consider two-loop correction 
to the propagator of the heavy scalar in original relativistic theory depicted in FIG. 3 a). The 
result of this diagram is proportional to the following integral: 

J^. = / : „ -„ . (9) 



[qf + ie] [ql + it] [{j) + qiY - + ie\ [{p + qi + ga) - + ie 
From the method of dimensional counting it follows that 



where 



k=0 k=0 k=0 



M2 - p2 -2vl 1 



p2 M 1 + ^ 



'11^ 



M 
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FIG. 3. Two-loop correction to the heavy scalar propagator in the heavy scalar-massless scalar 
scattering process. 



and the coefficients fik are determined by the original integral (|]) [T^. Substituting Eq. (pJ]) 
into Eq. (|I0[) we obtain: 



J,, = {-2vir-'M-'±D,,{'^-^] +{~2vir-'M--'±D, ^^^'^ 



A;=0 



^2k 



k=0 



M 



k=0 



M 



(12) 



where Dik do not depend on M, I or v. 

Actual calculations of J21 can be performed using the methods of as follows. Let us 
rewrite: 



J: 



21 



[?i + ^e] + ^e] [^Mv ■ (/ + gi) + (/ + qif + ie] [2Mv ■ {I + qi + q^) + {l + qi + ^2)^ + ie 

(13) 

First we expand the integrand in inverse powers of M, change the order of integration and 
summation and obtain: 



J. 



(1) 



21 



dP-qidr-q2 



4M3 J [ql + it] [ql + it] [v ■ {I + qi) + ie] [v ■ {I + qi + ga) + ie] 



+ 



4M4 \J [g2 + [^2 ^ [t; ■ (/ + gi) + ie]^ [v {I + qi + ga) + «e] 
1 r rqirq2 {qi + q2 + lf ^ +■■. 



2 J [ql + ze] [ql + ie] [t; ■ (/ + gi) + ie]' • (/ + gi + ga) + it] 



(14) 
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Second we re-scale q2 q2M, extract a non-integer power of mass, expand the integrand 
in inverse powers and change the order of integration and summation. The result is: 



t(2) 
^21 



M 



n— 6 



d"-qid"'q2 



4 J [g2 + ie] [ql + ie] [v ■ {l + qi) + ief [ql + 2vq2 + ie 



M 



n~7 



rqirq2 {qi + lf 



+ 



4 bl + ^e] bl + ^e] [v-{l + qi)+ ief [ql + 2vq2 + it] 

d^qid^q2 {v ■ {qi + I) - 2q2 •(?! + /)} \ ^ _ _ _ 

[q! + ie] [ql + ie] [v ■ {I + qi) + ief [ql + 2vq2 + ief J 



(15) 



and third we re-scale qi — > Mqi, q2 — > Mq2, extract non-integer power of the mass, change the 
order of integration and summation and obtain: 



7(3) 

•^21 



M 



2n-10 



d^qid"q2 



-M 



2n-ll 



[ql + ie] [ql + ie] [qf + 2v qi + ie] [{qi + ^2) + 2v ■ (gi + ^2) + ie 
d''qid''q2 {Al-qi+Avl) 



[ql + ie] [ql + ie] [ql + 2vqi+ ie] {qi + ga) + 2v ■ (gi + ga) + ie 



+ 



[ql + ie] [ql + ie] [ql + 2v qi + ie 



d"gid"g2 {21 ■ (gi + ga) + 2v ■ 1} 

' ■ ' h + q2f + 2v ■ {qi + q2) + ie 



+ ■ 



(16) 



Integral J21 is nothing else than a sum of J2i\ J21 and J21 

Evidently, J21 is expandable in inverse powers of M and hence heavy baryon approach 
reproduces this part of J21. ^21^ can not be reproduced (because it contains M^""^^) but it 

is analytic in / and hence it can be taken into account by adding compensating terms into 

f 2) 

the Lagrangian of the heavy baryon approach. It is J21 , corresponding to the second term in 
Eq. (p!2D which is non-trivial and might cause problems: there can appear terms which are not 
expandable in powers of M and have non- analytic dependence on /. 

This feature does not appear on a one loop level. Let us consider this problem in details. 

The above given representation for J21 can be obtained also as follows: One loop sub- 
integral over g2 can be represented as a sum of two parts: the first part is a result of expanding 
integrand of this sub-integral in inverse powers of M and changing the order of integration and 
summation. The second part is obtained by rescaling g2 — > g2M, extracting non-integer factor 
of M, expanding the integrand in powers of 1/M and changing the order of integration and 
summation: 



r(3) 



J'. 



21 



[ql + ie] 2Mv ■ (/ + gi) + (/ + qif + le 



{Fi(M,/ + gi) + M"-^F2(M,Z + gi)} (17) 



where Fi and F2 represent series in 1/M. As we concluded from one- loop analysis heavy baryon 
approach reproduces Fi at one loop level and M"-~^F2 is reproduced by adding compensating 
terms into the Lagrangian of the heavy baryon approach. 
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FIG. 4. Two-loop correction to the heavy scalar propagator in the heavy scalar-massless scalar 
scattering process. 



Now, expanding the denominator appearing in Eq. (17) 

1 



[qf + ie] 2Mv ■ (/ + qi) + (/ + qif + ie 



in inverse powers of M and changing the order of integration and summation in Eq. ([T7|) we get 
the result which is equal to + J21 . This makes clear that heavy baryon approach reproduces 
J21'* which was addressed as a possible source of the trouble. As for J21'' it is analytic in / and can 
be reproduced by compensating terms. Note that J21 is obtained from Eq. (|1^) by rescahng 
qi — i> qiM: one extracts non-integer power of M, expands the integrand in powers of 1/M and 
changes the order of integration and summation. Doing so one gets M"~^/i + M^"~^°/2, where 
M'^n-io is equal to J^i and /i turns out to be a sum of trivial terms (zeros). 

FIG. 3 b), FIG. 3 c) and FIG. 3 d) correspond to J21 ■, J^i and respectively. 

Next we consider two-loop correction to the heavy scalar propagator in original relativistic 
theory depicted in FIG. 4. The result of this diagram is proportional to the following integral: 



J, 



22 



J [ql + ie] [ql + ie] [{p + qif - + ie 



(i"gi(i"'g2 

{p + qi + q2f -M^ + ie 



{p + (12? -M^ + ie 



For later use before analysing Eq. ([T8|) let us consider Jy - an off-mass shell integral of the 
one loop correction to the vertex: 



d'^q 



[g2 + te] [{p + g)2 - + te] [{k' + qf - + ie 
where p = Mv + I and k' = Mv + I'. 

dJ'q 



(19) 



[2Mv ■{l + q) + {l + qf + ie] [2Mv ■ {V + q) + (/' + qf + ie] 



JUlJ') + J^ilJ') (20) 
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where is obtained by expanding the integrand of Eq. (pO]) in inverse powers of M and 
changing the order of integration and summation: 



4M2 J [q'^ + ie][v-l + v-q + ie\[v-l' + v-q + ie] 
Iff d^'q 



+ / ,]+■■■ (21) 

J [q"^ + ie] [v ■ I + V ■ q + ie] [v ■ I' + V ■ q + ie] J 

is obtained by rescahng q —>■ qM, extracting non-integer power of M, expanding the 
integrand in negative powers of the mass and changing the order of integration and summation: 

J^il, n = M"- / ; - - M"- d-q2{v.l + l.q) 



'" ' [g^ + ie] [2t> ■ g + g2 + ze]^ \J [g^ + ie] [2t> ■ g + + ie]^ 

(22) 



[g^ + %t\ [2t> ■ g + g2 + ie] 

Heavy baryon approach reproduces J^. is analytic in / and /' and can be reproduced by 
adding compensating terms into the Lagrangian of the heavy baryon approach. 



Applying the method of dimensional counting to J22 we obtain the following expression: 




A:=0 \ / fc=0 \ / 

(23) 



where A^^ do not depend on M, I or v. In Eq. ( [23| ) J22 is the result of expanding integrand 

(2) (3) 

in inverse powers of M and integrating the series, J22 and J22 are obtained by rescaling 
gi ^ giM and g2 — > g2M correspondingly with subsequent expansion of the integrand and 
change of the order of integration and summation and 722"* is the result of the simultaneous 
rescaling gi — qiM and g2 — > g2M and integration of the resulting expansion. 

Heavy baryon approach reproduces straightforwardly J22'*; J22 is analytic in momenta and 
hence can be reproduced by compensating terms in the Lagrangian of the heavy baryon ap- 
proach. The terms J22'' and J22'* are not expandable in 1/M and they are not analytic in 
momenta. In a full analogy with the previous analysis for J|j^ these terms are reproduced by 
taking into account the contributions of the compensating terms which have to be introduced 
into the Lagrangian of the heavy baryon approach in order to reproduce the expression for the 
one loop sub-diagrams of this two-loop diagram. 
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To see that this is actually the case let us represent J22 in the following way: 



22 



m + ie] [{p + qiY - M2 + ie\ J [g| + le] [{p + gi + 52) - + ie\ [{p + gs) - + ie 



[gf + ie] 2Mv ■ {I + qi) + (/ + qif + 



(24) 



Note that J^(/ + gi, /) corresponds to compensating terms included into the Lagrangian of the 
heavy baryon approach. 
Expanding denominator 



1 



2Mv ■ {I + qi) + (/ + qiY + ie 



(25) 



in 1/M and changing the order of integration and summation we reproduce 722"* + ^22''- So, 
heavy baryon approach reproduces these two terms (here J22^ occurred because we included 
contributions of compensating terms corresponding to one-loop sub-diagrams). 

J22^ and J22'' ^'I's reproduced by rescaling qi —>■ qiM, extracting non-integer factors of M, 
expanding integrand in 1/M and changing the order of integration and summation. 

J22^ is analytic in momenta and hence can be reproduced by compensating terms included 
into the Lagrangian of the heavy baryon approach. As for J22 it is equal to J22 and comes 
from one-loop compensating terms as well. This fact should be clear from FIG. 4 where FIG. 4 
b), FIG.4 c), FIG.4 d) and FIG.4 e) correspond to 4^, 4^, 4f and 4^ respectively. 
Analogous results are obtained for all the remaining two-loop diagrams. 
From the above analysis it follows that heavy baryon approach reproduces the results of the 
original relativistic theory at a two-loop order. 



r(2) 



IV. CONCLUSIONS 

In this work we have addressed the problem of matching of heavy baryon approach to 
the original relativistic theory. The heavy Baryon approach corresponds to the expansion 
of the integrand in inverse powers of the large mass with subsequent change of the order 
of integration and summation. As this two procedures are not commutative, the difference 
has to be compensated by adding terms into the Lagrangian of the heavy baryon approach. 
As the addressed problem does not actually depend on the details of the given model we 
considered a simple example of heavy and massless interacting scalar fields. Using the method 



of calculation of loop integrals by dimensional counting outlined in we analysed one and 
two loop diagrams and demonstrated how the difference between relativistic and heavy baryon 
calculations is compensated by adding terms to the Lagrangian of the heavy baryon approach. 
At two-loop level the difference can be compensated only after one includes the contributions 
of compensating terms for one loop sub-diagrams. While we included only selected diagrams 
in this paper, the very same conclusions are valid for one and two loop diagrams which were 
not included in here. 
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We believe that the iterative procedure of considering contributions of compensating terms 
for one-loop diagrams in two-loop calculations which is crucial to resolve the matching problem 
leads to analogous results for higher loops. 

While we considered a simple model of scalar fields the problems of interchange of integra- 
tion and expansion in inverse powers of heavy particle mass are the same for more realistic 
models with included fermionic and vector fields. In heavy baryon chiral perturbation theory 
the compensating terms with similar structure are actually summed up and included as redef- 
initions of already existing coupling constants. This redefinition is crucial, it actually leads 
to the consistent power counting of the heavy baryon chiral perturbation theory. In heavy 
baryon approach the coupling constants which correspond to re-defined relativistic coupling 
constants are introduced as initially free parameters which are to be fixed from experimental 
data. Working up to some given order in heavy baryon approach one actually re-sums low 
order contributions of an infinite number of relativistic high order loop diagrams. 

It was shown in Ref. that, within HBCHPT, an infinite number of internal line insertions 
must be summed to describe the scalar form-factor of the nucleon near threshold. As we 
demonstrated above the heavy baryon expansion reproduces the expansion of the relativistic 
result. This conclusion is formally still correct for the scalar form-factor of the nucleon, but 
the problem is that the expansion of the relativistic result is not convergent near threshold. 
This problem has been successfully resolved recently by Becher and Leutwyler using "infrared 



regularization" |12 



As was demonstrated above relativistic diagrams contain parts which can not be altered by 
adding local terms into the Lagrangian. These parts are directly reproduced by heavy baryon 
approach and they respect power counting. Other parts which are responsible for violation of 
the power counting in relativistic theory can be changed by adding counter-terms. Hence it 
should be more or less clear that the problems of the relativistic approach, in particular that 
multi-loop diagrams contribute into low order calculations encountered in [^, can be solved 
within relativistic approach using appropriately chosen normalisation condition. Hence one 
could from the very beginning work within original relativistic approach and never encounter 
the problems of near threshold behaviour of the scalar form-factor of the nucleon. These 
problems will be addressed in next paper. 
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